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Abstract: We develop numerical algorithms for solving the Einstein equation on Calabi-Yau 
manifolds at arbitrary values of their complex structure and Kahler parameters. We show 
that Kahler geometry can be exploited for significant gains in computational efficiency As a 
proof of principle, we apply our methods to a one-parameter family of K3 surfaces constructed 
as blow-ups of the T 4 /Z2 orbifold with many discrete symmetries. High-resolution metrics 
may be obtained on a time scale of days using a desktop computer. We compute various 
geometric and spectral quantities from our numerical metrics. Using similar resources we 
expect our methods to practically extend to Calabi-Yau three-folds with a high degree of 
discrete symmetry, although we expect the general three-fold to remain a challenge due to 
memory requirements. 
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1. Introduction 

In 1977, S.-T. Yau [jjj proved E. Calabi's conjecture |2| that a compact Kahler manifold with 
vanishing first Chern class admits a Ricci-flat metric in each Kahler class. Interest in Calabi- 
Yau manifolds was subsequently generated among physicists by the discovery that they can 
serve as supersymmetry-preserving compactification manifolds in string theory ||. In the 
three decades since the proof of Yau's theorem, many examples of Calabi-Yau manifolds have 
been constructed and studied, principally using methods of algebraic geometry, and much has 
been learned about their mathematical properties and physical applications. Nonetheless, a 
major gap in our knowledge about Calabi-Yau manifolds has persisted during this time, 
namely the Ricci-flat metrics themselves. Yau's proof is not constructive, and no example of 
a smooth Ricci-flat metric is known explicitly for any Calabi-Yau manifold. Indeed, perhaps 
we should not expect there to exist such metrics in closed form 0]. 
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The question thus arises as to whether it is possible to solve the Einstein equation nu- 
merically on a Calabi-Yau manifold. The purpose of this paper it to show that it is. We 
describe a method for doing so, and display the metrics obtained by applying that method 
to the smallest-dimensional Calabi-Yau, the K3 surface. 

The algorithms we developed rely in an essential way on the underlying complex and 
Kahler geometry of Calabi-Yau manifolds. In fact, one of the main points of this paper is to 
show that those properties are as powerful for numerical work as they have already proven 
to be for analytical calculations and proving theorems. To explain this, let us consider the 
challenges faced by someone attempting to numerically solve the Euclidean Einstein equation 
on a general real four-manifold. In some sense this is a problem in numerical relativity, and 
one can get a sense of its scale by considering that four-dimensional problems in numerical 
relativity, such as black hole collisions, can be solved (if at all) only with the investment 
of extremely large computing resources, typically supercomputers. (Note that Calabi-Yau 
manifolds admit no continuous isometries that could reduce the effective dimensionality of 
the problem.) At a more fundamental level, whereas conventional numerical relativity deals 
with solving the Einstein equation as an initial-value problem on a Lorentzian spacetime, 
here we wish to solve it on a Euclidean manifold, a problem for which general algorithms are 
lacking. The challenges for creating an algorithm include the usual issues of gauge fixing and 
coordinate singularities, as well as finding a way to fix any moduli the solutions might have. 
An even more difficult challenge would be to avoid (or deal with) the curvature singularities 
that generically form under relaxation schemes such as Ricci flow. 

Let us now see how the framework of complex and Kahler differential geometry allows 
one to naturally solve, or greatly ameliorate, each of these problems in turn. Firstly, the 
Kahler formulation of geometry can be employed to vastly reduce the scale of the problem 
compared to the language of real differential geometry. The metric can be encoded in a 
single scalar function, the Kahler potential (as reviewed in subsection 2.1). In terms of the 
Kahler potential the Einstein equation takes the relatively simple form of a Monge-Ampere 
equation (as reviewed in subsection 2.2). So we have a large simplication of the many degrees 
of freedom required to locally parameterize a general metric, and of the complicated set of 
differential equations in the usual form of the Einstein equation. 

Secondly, complex coordinates offer a naturally adapted gauge choice, since with respect 
to any Kahler metric they satisfy the harmonic gauge condition, a well-known gauge choice for 
numerical relativity. Furthermore the gauge fixing is nearly complete: there are only a finite 
number of (continuous) residual gauge transformations, since compact complex manifolds 
admit only a finite number of holomorphic vector fields. In fact, for Calabi-Yaus this number is 
zero, so the gauge fixing is complete. Since there are no gauge transformations, no coordinate 
singularities can appear. 

Thirdly, the moduli of Calabi-Yau manifolds, which are divided into complex structure 
and Kahler moduli, may be fixed at any desired values — before solving the Einstein equation — 
in the following way. With the manifold defined topologically by an atlas of patches, the 
complex structure is fixed by fixing the holomorphic coordinate transition functions on the 
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patch overlaps, while the Kahler moduli are fixed by fixing the Kahler transformations on 
them (which then serve as the boundary conditions for the Kahler potential). (Details of the 
above are given in subsection |2.1| .) 

Finally, there is a kind of stability that appears to be inherent in Kahler geometry, 
miraculously eliminating the problem of spontaneous formation of curvature singularities. 
Ricci flow, for example, contrary to its behaviour on real manifolds, is extremely robust 
on Kahler manifolds, as shown by Cao's long-time existence theorems ||. On Calabi-Yau 
manifolds in particular, starting from any Kahler metric it converges to the Ricci-flat metric 
in the same class. In principle therefore it provides a general algorithm for solving the Einstein 
equation on Calabi-Yaus. However, since it is a rather inefficient method from a computational 
viewpoint (rather like solving the Laplace equation by simulating diffusion), we developed and 
used instead a Gauss-Seidel-type relaxation algorithm for the Monge- Ampere equation. While 
we lack a convergence theorem for our algorithm, we found that in practice it was just as 
robust as Ricci flow, presumably for the same underlying reasons. (In fact, in one sense our 
algorithm was even more robust than Ricci flow, since remarkably it converged even when 
the initial Kahler potential didn't define a [positive-definite] metric.) Subsection 2.3 contains 
a discussion of Ricci flow and a description of our algorithm. 

As a proof of principle for our methods, we applied them to a class of K3's known as 
Kummer surfaces, which are blow-ups of the orbifold T 4 /Z2- For simplicity, we considered 
the most symmetrical Kummer surfaces, namely those for which the torus is cubical and all 16 
fixed points are blown up identically. This symmetry leaves only one modulus (not counting 
the trivial volume modulus), namely the ratio of the size of the blow-ups to the size of the T 4 . 
In fact there is a limit to how large this ratio may be, since at some point certain holomorphic 
curves shrink to zero size, signalling the appearance of new orbifold singularities — on this 
wall of the Kahler cone the manifold is an orbifold of another, smooth K3. (Details of the 
construction are given in subsection [O].) Using our algorithm, we computed the Ricci-flat 
metric at various points over the full range of this modulus. Each point required a few days 
on a garden-variety desktop computer for the highest resolutions. Subsection |3.2| is devoted 
to an exploration of the resulting geometries as a function of the modulus. Various curvature 
invariants are plotted, as well as a low-lying eigenvalue of the scalar Laplacian. 

The success we had with these highly symmetrical K3 surfaces leads us to consider pos- 
sible generalizations. Using our methods, could we solve the Einstein equation on generic 
K3's, which lack such discrete symmetries, or on Calabi-Yau three-folds? What about on 
other Kahler manifolds such as del Pezzo surfaces, or on Calabi-Yaus with matter such as 
fluxes and branes? The estimates we make in subsection 4.1 suggest that the answer to these 
questions is yes, but that due to memory limitations in the three-fold case we would need to 
be helped by a high degree of discrete symmetry. 

In subsection 4.2 we return to the problem of solving the Einstein equation in the real 
Euclidean context, and explore what this work has taught us that might generalize to non- 
Kahler geometries. We conclude in subsection |4.3j by mentioning some possible mathematical 
and physical applications of these numerical metrics. 
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The C code for our simulations, as well as animated versions of the plots shown in this 
paper, are available at the website http://schwinger.harvard.edu/~wiseman/K3/, 



2. General method 

The purpose of this section is to discuss in a general way the problem of solving the Einstein 
equation numerically on a Calabi-Yau manifold at a given point in its moduli space. We first 
briefly review the essentials of Kahler geometry, 1 describing how the geometry is encoded in 
the Kahler potential and how the moduli are fixed. We then explain how, in Kahler geometry, 
the Einstein equation reduces to a Monge-Ampere equation. Finally, we discuss in general 
terms the numerical algorithms we applied to solving that equation. 

2.1 Kahler geometry 

We work on a manifold with a fixed complex structure, that is, on each coordinate patch U a 
we have a set of complex coordinates {z^z^} such that on each overlap U a n Up the transi- 
tion functions Zp(z a ) are holomorphic (a, (3 index the patch and i,j the complex coordinate). 
Fixing the complex structure restricts the allowed gauge transformations to holomorphic dif- 
feomorphisms. Continuous holomorphic diffeomorphisms are generated by holomorphic vector 
fields. In contrast to the infinite number of vector fields generating general diffeomorphisms, 
it can be shown that a compact manifold admits only a finite number of holomorphic vector 
fields. Even better, Calabi-Yau manifolds have none at all, for the following reason. Every 
Calabi-Yau is equipped with a nowhere-vanishing holomorphic (n, 0) form Q (which depends 
on the complex structure but not on the metric). If v % is non- vanishing with holomorphic 
components, then the same is true of the (n — 1,0) form v ll £li 1 i 2 ...i n dz %2 A • • • A dz ln . Since 
the Hodge number /i n_1, ° vanishes on a Calabi-Yau, such a form must be d exact, which is 
impossible for an (n — 1, 0) form. Hence, on a Calabi-Yau, fixing the transition functions for 
the complex coordinates amounts to a full gauge fixing. As mentioned in the introduction, 
this is a great advantage for numerical work. 

A metric is called Kahler with respect to a given complex structure if it is Hermitian, 
9ij = 9i] = 0, and if the associated Kahler form J = ig%jdz % A dz 3 is closed, dJ = 0. J 
is obviously a positive (1, 1) form, since the metric is everywhere positive definite. The 
cohomology class of J is called the Kahler class. The set of potential Kahler classes, i.e. (1, 1) 
cohomology classes containing at least one positive form, is called the Kahler cone, or Kahler 
moduli space. 

The fact that the Kahler form is closed implies that it, and the metric, are locally ex- 
pressible as the matrix of second derivatives of a scalar, 

9ij\u a = d i d 3 K °" J \u a = iddKa, (2.1) 

A detailed review of Kahler geometry and Calabi-Yau manifolds may be found in the excellent set of 
lecture notes [[|. 
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where K a is a real function, the Kahler potential, defined on the patch U a . Given gij, the 
Kahler potential is not unique but can be changed by adding the real part of any holomorphic 
function, a so-called Kahler transformation. From a numerical point of view, one advantage 
of using the Kahler potential to encode the geometry is obvious: it reduces the number of 
functions to store from n(2n + 1) for the full metric (or n 2 if the hermiticity condition is 
imposed) to a single one. 

The volume of the manifold is given in terms of the Kahler form by 

V=— t I J n . (2.2) 
n\ J 

Therefore on a compact manifold J cannot be exact, or else the manifold would have zero 
volume. So the Kahler potentials obtained by integrating (^[l]) must disagree on the patch 
overlaps by Kahler transformations: 

K a - Kp = u a/3 . (2.3) 

These Kahler transformations serve as boundary conditions for the Kahler potential (the only 
boundary conditions if the manifold is compact). In doing so, they perform two important 
tasks. First, they fix the Kahler class. To see this, note that two Kahler potentials {K a } and 
{K' a } which have the same Kahler transformations, K a — Kp = K' a — K'a, differ by a globally 
defined real function <j). Therefore the difference between the corresponding Kahler forms, 
J — J' = iddcf), is an exact form. Conversely, given a representative of some particular Kahler 
class, a corresponding set of u a ps may easily be found by solving ( ^l| ) separately on each 
patch. (The representative need not be positive. This is useful particularly near the edge of 
the Kahler cone, where it may not be easy to find positive representatives.) The second task 
performed by the Kahler transformations is to (almost) eliminate gauge-equivalent Kahler 
potentials, that is, different Kahler potentials that give rise to the same metric: for each 
Kahler form J in the class defined by a set {u a p}, there is a unique solution to fl2.1|) and 
( [2. 3D , up to constant shifts of all the K a . 

When the metric is Kahler, only the purely holomorphic and antiholomorphic components 
T l - k and T l _^ of the Christoffel symbol are non-zero. Together with the hermiticity of the metric, 
it follows that the coordinates z 1 and z l are harmonic, a well-known gauge fixing condition for 
numerical relativity. However, Kahlerity is a stronger condition than harmonicity; whereas a 
given metric always admits a harmonic coordinate system, at least locally, a generic metric 
is not Kahler with respect to any coordinates, even locally. 2 

To summarize: We fix the complex structure by specifying the coordinate transition 

functions on the patch overlaps. Using the Kahler potential to encode the geometry, we fix 

the Kahler class by specifying the Kahler transformations on the overlaps. With this set-up 

in hand, we now turn to the question of what equation to solve for the Kahler potential. 

2 Nor are Kahler coordinates necessarily unique for a given metric. For example, a Ricci-flat metric on a 
hyperkahler manifold, such as K3, is Kahler with respect to a continuous family of different complex structures. 
We will not make use of this additional structure, but will be content to fix a particular complex structure at 
the outset. 
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2.2 The Monge- Ampere equation 

There exists a simple expression for the Ricci tensor of a Kahler metric: 

R hl = % = 0, R kJ = R Jk = -d k dj\n det g i3 (2.4) 

(note that det ^ = Just as we defined the Kahler form J, we can define the Ricci 

form 1Z = iR k fdz k A dz l . By virtue of ( |2.4[) , 1Z is closed (but not necessarily exact, since 
det gij is not in general a globally defined function). One can show that its cohomology class, 
the first Chern class ci, is a topological invariant: two different Kahler metrics will give rise 
to Ricci forms that differ by an exact (1,1) form. Calabi-Yau manifolds are defined by the 
condition that c\ vanishes, in other words that the Ricci form is always exact. According to 
Yau's theorem jlj], on a Calabi-Yau manifold each Kahler class (i.e. (1,1) cohomology class 
containing at least one Kahler form) contains a unique Ricci-flat Kahler form. 
Prom fl2.4p , the Einstein equation for a Kahler metric reads 

dfcC^lndet = 0. (2.5) 

In terms of the metric this is a second-order PDE, but in terms of the Kahler potential 
it is fourth order. It might seem then that working with the Kahler potential is not so 
advantageous after all. However, as we will now explain, by the magic of complex analysis 
it can be reduced to a second-order PDE, specifically a Monge-Ampere equation (a PDE in 
which the derivatives appear in the form of a Hessian). 

For simplicity, let us assume first that the coordinates have been arranged in such a way 
that the Jacobians detij^dz^/dzp) of the transition functions are 1 on all the overlaps. In that 
case det is a globally defined function, and on a compact manifold the Einstein equation 
([2.5D is equivalent to it being constant: 

det = A. (2-6) 

This is a non-linear Monge-Ampere equation for the Kahler potential. The constant A is 
related to the volume of the manifold, given by Q2.2D , which depends only on the Kahler class 
and may therefore be calculated a priori from the Kahler transformations on the overlaps (see 
Appendix [B| for details). 

The coordinate system used for the Kummer surface in Section ^ happens to satisfy the 
condition of unit coordinate Jacobians assumed in the previous paragraph. It is instructive 
nonetheless to consider what one would do in the more general situation. Clearly we need to 



replace the right-hand side of (2J3) with something that has the correct transformation law 
on the overlaps, and also that implies the Einstein equation (|2.5| ) . The holomorphic (n, 0) 
form Q once again comes to the rescue; the generalization we seek is: 



det0ij = A|fii.„„r (2.7) 

(this can also be written in the coordinate-invariant form J n = (—i) n n\ AO AH). If Q is not 
known explicitly, then |r2i... n | 2 can be found as follows in terms of an arbitrary Kahler metric 
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(jij (which need not be in the same Kahler class as the desired solution, since depends 
only on the complex structure). Since the manifold is Calabi-Yau, its Ricci form is exact, 
Rij = didjF, where F is a globally defined function which can be calculated explicitly (either 
analytically or by standard numerical methods). We then have 

|^i...n| 2 = e p det gij, (2.8) 

and equation (^?j) may be applied. 

To get a feeling for the Monge- Ampere equation, it is useful to consider equations Q2.7|) 
and (|2.8[) in the case where gij is in the desired Kahler class. Then we can write 

K a = K a + <j), (2.9) 

where <j> is a globally defined scalar. In terms of <j) the Monge- Ampere equation is 

det($ + g k ~ 3 didj(f)) = \e p . (2.10) 

The left-hand side is a non-linear operator acting on <f>. If didj4> is small (i.e. if g^ is almost 
Ricci-flat) then we can linearize, yielding a Poisson equation: 

^V 2 + O {didjcj)) 2 = \e p - 1. (2.11) 

In this sense Yau's theorem can be understood as a generalization of the existence theorem 
for the Poisson equation. 

2.3 Methods 

The similarity to the Poisson equation observed in the last subsection suggests that methods 
for solving it might be generalized to the Monge-Ampere equation. Among the standard 
methods for the Poisson equation, some are local and others are non-local. We have restricted 
ourselves to local schemes, for two reasons. First, the Monge-Ampere equation is local and 
non- linear. Second, for most manifolds (even highly symmetrical ones), the individual patches 
U a (or, more precisely, their images in C n , i.e. the ranges of the coordinates z^), have rather 
irregular shapes. This motivated the use of a lattice discretization, rather than non-local 
spectral representations such as Fourier modes or wavelets. 

One simple (but inefficient) method for solving the Poisson equation is to simulate dif- 
fusion. The analogous geometric relaxation equation is Ricci flow, which is defined by a first 
order equation in an auxiliary time dimension: gij = —R^. This flow is well studied mathe- 
matically H, largely because of its application to geometrization. It is also of some interest 
in physics as the one-loop renormalization-group evolution for the target space geometry of 
a sigma model. In terms of the Kahler potential the flow is governed by 



Note that on the overlaps K a — Kr = 0, so the Kahler transformations, and hence the Kahler 
moduli, are conserved. Cao || has shown that Ricci flow starting from an arbitrary Kahler 
metric converges to the Ricci-flat metric in its class. 

Some work has been done on numerical simulation of Ricci flow on real geometries 
|lC| , [TT| , but not (as far as we know) in the Kahler case. We experimented with simple lattice 
implementations of such a scheme. However, if one is only interested in the endpoint of 
the flow, namely the Ricci-flat metric, as opposed to the whole flow, this method is clearly 
very inefficient, as it requires solving an equation in one higher dimension. Furthermore, 
stability of the diffusion problem typically requires implicit finite differencing schemes which 
are rather inconvenient (particularly in several dimensions). Another more subtle drawback 
which nonetheless is rather serious in practice is that it requires an initial Kahler (i.e. positive) 
form in the desired class, as seen above explicitly from the logarithm in ( |2.12| ). As mentioned 
in subsection 2.1 , it is not always easy to find a positive representative of a given class, 
especially near the edge of the Kahler cone. 

The prototype local method for solving the Poisson equation is the Gauss-Seidel method. 
Here the discretized Poisson equation is solved at each lattice point in turn. After a suitable 
number of iterations over the whole lattice, the discretized equations will be solved to a given 
accuracy. This is very robust and simple to implement. Furthermore whilst Gauss-Seidel is 
slow compared to spectral or multigrid methods in low dimension, scaling as N l+2 / d rather 
than iVlog N in d real dimensions (where N is the number of lattice points), in our K3 case 
of 4 real dimensions, and more so for still larger dimensions, the advantage is not so great. Of 
course multigrid could be implemented relatively simply to improve our Gauss-Seidel method 
if speed became a crucial issue. 

Our analog of the Gauss-Seidel method for the Monge- Ampere equation is as follows. On 
a lattice the metric is determined from the Kahler potential by taking discrete derivatives. 
Thus the value of det at a given site is a function of the values of K a at that site and 
its neighbors, out to some distance depending on the order of finite differencing used. Our 
algorithm directs one to go through each site of the lattice in turn, changing K a at that site 
to the value which solves (|2.7|) given its values at the neighboring sites. Although we do 
not have a convergence theorem for this algorithm, we found that in practice it converged 
on the full range of Kahler parameters studied. This included cases where the initial Kahler 
potential did not define a positive (1, 1) form. 

As noted earlier, the constant A which relates the coordinate volume to the proper volume 
can be computed analytically from the Kahler class. Now consider solving the equation 

det g i3 = A. (2.13) 

We will only find a solution to this new Monge- Ampere equation when we set A = A. However, 
there is a subtlety, as when we implement this equation numerically, the errors involved in 
finite differencing will mean that the value of A that solves the finite differenced equation will 
actually differ slightly from the true continuum value A. This implies the finite differenced 
equation should have no solution when A = A, which sounds mildly disastrous. However, 
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we can see in detail how this discretization error manifests itself if we consider the Ricci 
flow equation ( 2.12|) with A replaced by A. Then, since the Monge-Ampere equation has no 



solution for A 7^ A, the flow will never reach a fixed point. However, the flow does asymptote 
to one with a simple time dependence, namely, 

K~ x (t) = K x + tv, v = \nt, (2.14) 

A 

where K\ is the solution of the Monge-Ampere equation for the true value of A. The constant 
v, as determined above, then gives the asymptotic time dependence of the flow. By analogy 
with the Ricci flow above, this implies that under Gauss-Seidel iteration, our finite differenced 
Monge-Ampere equation ( |2.13| ) will have an asymptotic solution that drifts in iteration time 
by a constant mode. This simply corresponds to a drifting Kahler transformation; therefore 
the real metric does indeed tend to a static solution, and disaster is avoided. However, since 
we would rather have a fixed endpoint to our Gauss-Seidel method, we "cure" this drifting 
due to discretization error by solving the Monge-Ampere equation fl2,13|) , and dynamically 
determine A by averaging det g^ as we perform the Gauss-Seidel iterations. This procedure 
then yields a static end solution, and the value A should approach the true analytic A as the 
continuum is approached by increasing the lattice resolution. 

It is important to check that the solution to the lattice version of the Monge-Ampere 
equation converges to the continuum solution as the lattice resolution is improved. For 
this purpose it is useful to have some quantities that can be calculated from the lattice 
solution, and for which the exact, continuum value is also known. Here we will mention 
three such quantities. The first is that mentioned directly above, namely the total volume 
of the manifold, in other words how close A is to A. The agreement between the analytic 
and numerical values tests not only the quality of the solution, but also the error in fixing 
the Kahler class. A second quantity that is known exactly is the Euler number x °f the 
manifold (which of course does not depend on the Kahler class). This can be calculated 
numerically in terms of the Kahler potential via a Gauss-Bonnet theorem. Finally, there is 
the difference between Kahler potentials K a — Ka on the overlap U a fl Up. In subsection 



2.1 it was argued that this difference should be set equal to the fixed Kahler transformation 
as a boundary condition for the K^s. More precisely, however, the boundary condition 
is imposed only on the edge of each patch. For the continuum equation, this is enough to 
guarantee K a — Kp = u a p also in the interior of U a nUp, by the uniqueness of the solution to 
the Monge-Ampere equation. The lattice will introduce an error into this equality; conversely, 
that error measures how well the lattice solution approximates the continuum one. 

3. Application to Kummer surfaces 

For a first application of the techniques described in the previous section, we turned to the 
lowest-dimensional Calabi-Yau manifold, the K3 surface. K3 has played important roles in 



algebraic and differential geometry, as well as in string theory. (The lecture notes [12] provide 
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an excellent review of the role of K3 in string theory. See also the more mathematically- 
oriented review p3| .) The moduli space of Ricci-flat metrics on K3 is 58-dimensional: 40 
complex structure moduli and 20 Kahler moduli, minus 2 for the hyperkahler identifications. 
One of the moduli is the overall volume of the manifold; the other 57 are non-trivial in the 
sense that they do not act on the Ricci-flat metric by any straightforward transformation. We 
studied a particular one-parameter family of K3's which admit a simple construction and have 
a high degree of discrete symmetry, which serves to reduce the number of lattice points to be 
simulated at any given resolution. In the first subsection below, we explain the construction 
of these K3's. In the next subsection we describe the results obtained in our simulations, 
giving several examples of the kind of concrete geometrical information that is available from 
the explicit form of the metric. 

3.1 Construction 

Among the simplest K3 surfaces to describe explicitly are the so-called Kummer surfaces, 
which are the orbifold T 4 /Z2 with its 16 singular points blown up. After explaining this 
construction, we will specialize to the ones with the largest discrete symmetry group, namely 
those constructed from a cubical T with all 16 singular points blown up to the same size. 
The only free parameters are the size of the T and the size of the blow-up. We will see that 
the Kahler cone defines a finite range of values for their ratio. 

As a complex manifold, the torus T 4 is parametrized by a pair of complex coordinates 
{z l ,z 2 ) G C 2 , identified under translations by a set of four linearly independent vectors 
("I,OeC 2 (a = l,...,4), 

(*V)~(*V) + (t>2,«2). (3.1) 

Choosing the vectors (v\,v 2 ) fixes the complex structure of the torus (there are equivalences; 
the complex structure moduli space is actually only 8 real dimensional). The parity map 

(z\z 2 )^(-z\-z 2 ) (3.2) 



is compatible with the equivalences (3.1), so we can quotient the torus by it. Furthermore, 



it acts holomorphically, so the result, T /Z2, is a complex manifold. More correctly, it's an 
orbifold, since the parity map has fixed points. There are 16 of them, located at (z ,z 2 ) = 
1 ^a na ( v ai v a) with n a = or 1, and each carries an A\ (or C 2 /Z2) type singularity. 

We obtain a smooth complex manifold, known as a Kummer surface, by blowing up the 
fixed points. Consider for example the fixed point at the origin. To blow it up, remove 
that point from the manifold and add two new patches, with coordinates (y, w) and (y', w') 
respectively, and transition functions (which are of course holomorphic) 

1 J.J%\ ( Z 1 I ~2\2 



M = \-,^'y ) = {^2 iz r)' (3 - 3) 

(z\z 2 ) = ±v^(y,l) = ±v^7(l,i/). (3.5) 
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To avoid complications, the ranges of the new coordinates should be bounded in such a way 
that they do not include the other fixed points. Each of the 16 fixed points of the orbifold is 
given its own (y,w) and (y',w') patches, with the same transition functions (3.3-|3.5|) except 
that (z 1 ,^ 2 ) is replaced by {z l } z 2 ) — \ J2 a na ( v ai v a)- Hence we have a total of 33 patches. 
There are three important points to note about the new (y, w) and (y ; , w') coordinate systems. 
First, the identification under the orbifold action ( |3.2| ) is automatic in them. Second, the 
origin has been replaced by the surface w = w' = 0, parametrized by y = 1/y'. This is a 
CP 1 (or S 2 ), and is homologically non-trivial; it is called the exceptional divisor. Finally, the 
transition functions ( |3.3| - |3.5|) all have unit Jacobian. Hence it is quite easy to write down the 
holomorphic (2, 0) form: 

n = dz 1 A dz 2 = dy Adw = dw' A dy' . (3.6) 

Kummer surfaces have 8 complex structure moduli (inherited from the T 4 ) and 20 Kahler 
moduli (the 4 inherited from the T 4 , plus the size of each of the 16 exceptional divisors). It 
can be shown that they are special cases of K3 surfaces. (The missing 32 complex structure 
moduli are due to the fact that we blew up, rather than deformed, the orbifold fixed points.) 

Both for simplicity and in order to reduce the number of lattice points simulated, it was 
advantageous for us to choose highly symmetrical Kummer surfaces. The T 4 was taken to be 
cubical; in other words, the periodicities ( |3,l| ) were given by 

1 - 1 + 1 ~ z 1 + i, z 2 ~ z 2 + 1 ~ z 2 + i, (3.7) 



z ~ z 



while the Kahler transformations are those obtained for the flat metric on a cubical T 4 of 
side length b: 



(3. 



(3.9) 



K{z 1 + l,z 2 )-K{z v ,z 2 ) = b 2 fw + i) , 

Kiz 1 +i,z 2 )-K(z\z 2 ) = b 2 (lmz l + ^\ , 

K(z\z 2 + 1) - K(z\z 2 ) = b 2 (Rez 2 + ^j , (3.10) 

K(z\ z 2 +%)- K(z\z 2 ) = b 2 (imz 2 + . (3.11) 

The coefficients on the four right-hand sides correspond to the four Kahler parameters of T 4 ; 
here since the T 4 is cubical, they are set equal to a common constant b 2 . Each blown up fixed 
point has only one Kahler modulus; without loss of generality the Kahler transformations 
may be taken as follows: 

(3.12) 
(3.13) 
(3.14) 





K(y,w) 


= a 2 


ln\z 2 


K {z^,z 2 ) - 




= a 2 


\ii\z 1 




K(y',w') 


= a 2 


ln|y| 
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All 16 fixed points were blown up to the same value of the modulus a 2 . 

How much discrete symmetry do these Kummer surfaces possess? Our choice of complex 
structure admits a holomorphic diffeomorphism group of order 2 s , all of which is respected 
by our choice of Kahler class. The generators of this group are as follows: the translations by 
the vectors (|, 0), (|, 0), (0, |), and (0, |), which generate a Z| group that maps the origin to 
each of the other fixed points; the rotation z 1 t— > iz , which generates a Z4 group; the diagonal 
rotation (z l ,z 2 ) (iz 1 ,iz 2 ), which (in view of the identification under parity) generates a 
Z2 group; and the exchange {z , z 2 ) 1— > (z 2 , z 1 ), which generates another Z2. Of course, these 
generators don't commute with each other, so the full holomorphic diffeomorphism group 
is complicated. In addition, the anti-holomorphic diffeomorphism (z 1 ,^ 2 ) ^ (z 1 ,^ 2 ) is a 
symmetry of the real metric. Uniqueness of the solution to the Monge- Ampere equation 
guarantees that every symmetry of the Kahler class is an isometry of the Ricci-flat metric in 
that class. Therefore it is sufficient for us to simulate only the fundamental domain of the 
symmetry group, a reduction potentially by a factor of 2 9 . 

As explained in Appendix [B[ the volume is easily calculated from the intersection matrix 
of the second de Rham cohomology group H2 (R) : 

V = l6 4 -4vr 2 a 4 . (3.15) 



As discussed in subsection [2^, the volume calculated numerically may be compared against 
this analytic result to give a measure of how accurately we are fixing our Kahler class. More 
details are given in Appendix [A]. 

Equation ( |3.15 ) shows that, as the size of the blow-up increases, it eats away at the 
volume of the manifold. This indicates an upper bound a 2 /b 2 < (8-7T 2 ) -1 / 2 . In fact, the real 
bound is slightly lower: a 2 /b 2 < (Air)" 1 . We first noticed this as an empirical fact in our 
numerical trials, but the reason is not hard to understand. The volume of a holomorphic 
submanifold depends only on the Kahler class, not the metric (the total volume being a 
special case of this). Therefore, a necessary condition for being inside the Kahler cone is that 
all the holomorphic submanifolds have positive area. Our symmetric Kummer surfaces have 
three types of holomorphic curves. The first type are the 16 exceptional divisors. From the 
Kahler transformation (|3. 14 ) restricted to the w = surface, one may show that each has 
area A = ira 2 . We thus have the condition a 2 > 0. Another holomorphic submanifold is the 
curve {z 1 = C} U {z 1 = — C}, which represents the same two-cycle for all values of C other 
than the special values C = 0, |, |, 5 + §■ Using (|3.1C , 3~l~l] ), its area is b 2 , so we have b 2 > 0. 
In the special cases C = 0,i,|,i + i the curve passes through points which are outside of 
the z patch, so one has to be more careful. As we show in Appendix [B], these represent 4 
different two-cycles; each has area 

A = h) 2 - 2vra 2 , (3.16) 

accounting for the above upper limit on a 2 /b 2 . Obviously, the same results apply to the curves 
of constant z 2 , so there are a total of 8 curves of this type. In fact they are rational curves, 



- 12 - 



since they have topology T 2 /Zi2 ~ S 2 . Each of them intersects each exceptional divisor at a 
single point, so what is happening as they shrink to zero size is that the exceptional divisors 
are so large that they touch each other and therefore can't be blown up any larger. 

An isolated rational curve shrinking to zero size implies the formation of an A\ orbifold 
singularity. In the next subsection we will explicitly confirm the formation of these orbifold 
singularities from our numerical solutions. The simultaneous formation of 8 A\ singularities 
naturally leads to the idea that the manifold in this limit is globally a Z2 orbifold. By counting 
the Euler number, one finds that it must be an orbifold of another, smooth K3. Indeed, it 



can be shown [jlj, 15] that it is a rather well-known K3, namely the Fermat quartic surface 



x 



f + x\ + x| + x\ = in CP 3 , in the Kahler class induced from the Fubini-Study metric 
on CP 3 . 3 Hence our one-parameter moduli space interpolates between the orbifold of T 4 at 
a 2 /b 2 = and the orbifold of the Fermat quartic at a 2 /b 2 = (47r) _1 . For our fixed complex 
structure, these endpoints of our modulus a 2 /b 2 represent the edges of the Kahler cone. 

It is worth remarking that an approximate analytical solution to the Einstein equation 
can be constructed in the limit a 2 jb 2 <C 1 by smoothly joining the Eguchi-Hanson metric (a 
Ricci-flat and asymptotically fiat metric on the blow-up of C 2 /Z2) onto the flat torus metric 
18|. As we will see in the next subsection, the numerical method does not perform as 



well in this regime of parameters because the manifold contains a region of high curvature, 
namely the vicinity of the exceptional divisors. Thus the numerical and analytic approaches 
are effective in complementary regimes. 

3.2 Results 

We applied the methods described in Section ^ to find the Ricci-flat metrics on the symmet- 
rical Kummer surfaces constructed in the previous subsection, at 9 different values of the 
modulus a 2 /b 2 . In terms of the combination 

a 2 

a = 4vr— , (3.17) 
tr 

which ranges from to 1, the points for which the metrics were calculated were 

a = 0.03, 0.13, 0.28, 0.50, 0.61, 0.72, 0.79, 0.85, 0.92. (3.18) 

(In figures [l], @, and ||] we have left out a = 0.79 for typesetting elegance.) Without loss 
of generality, the volume of the manifold was fixed to be 1. The computational aspects of 



3 It was shown in Jl4|] that (the blow-up of) the orbifold of the Fermat surface by the Z2 action 
(asi, X2, £3, £4) 1 — > (xi, X2, — £3, — £4) has the same complex structure as the Kummer surfaces we consider. 
This action has 8 fixed points, and the resulting 8 exceptional divisors are identified with our "shrinking ratio- 
nal curves". It can furthermore be shown that the Kahler class on the Kummer surface at a 2 jb 2 — (47r) _1 lifts 
to the Kahler class on the Fermat surface induced from the Fubini-Study metric on CP 3 E3|. Amusingly, it 
has also been shown Jl3| |l6| , that the sigma model with this Kummer surface (or equivalently the orbifold 
of the Fermat surface) as its target space, at a particular volume and equipped with a particular B-field, is 
dual to one whose target space is the Kummer surface at the opposite edge of the Kahler cone, a /b 2 = 0! We 
are grateful to K. Wendland for very helpful discussions on these issues. 
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the problem — lattice discretization, convergence, etc. — are discussed in detail in Appendix [A|. 
Let us simply mention here that, at each of the above values of a, the metric was computed 
at four different lattice resolutions, labelled A, B, C, D; B has twice the linear resolution (or 
16 times the total number of points) of A, and so on. 

In this subsection, we will illustrate the kind of concrete geometric information that is 
available once the Ricci-fiat metric is known. We will focus on three ways to characterize 
the geometry: the distribution of the Euler density; the induced geometry on the exceptional 
divisors and other rational curves; and the low-lying spectrum of the Laplacian. If unspecified, 
results presented are generated from the highest resolution metric, D. 

On a Ricci-flat manifold, the simplest non-trivial curvature invariant one can construct 
is the square of the Riemann tensor (sometimes referred to as the Kretschmann invariant). 
In four dimensions, this happens to be proportional to the Euler density p: 



X 



V9P, p = ^R iJkT R» kI . (3.19) 



Thus its integral over the manifold is fixed, in this case at 24. Figures |l] and ^ show surfaces 
of constant p on the three-dimensional slice Imz 2 = at eight different values of a. 4 At the 
smallest value (a = 0.03), the curvature is highly concentrated near the fixed point, and is 
spherically distributed; here the metric closely approximates the Eguchi-Hanson metric (for 
which the isosurfaces of p are spherical in these coordinates, although the full geometry is only 
axisymmetric) . As we increase a, the curvature spreads out from the fixed point. However, it 
does not diffuse uniformly over the manifold. Instead, starting in the a = 0.72 figure (third to 
last), it gathers along the z 1 = and z 2 = planes. These are the rational curves, discussed 
in the last subsection (and referred to as c\ and C2 in Appendix |B|), that shrink to zero size 
as a approaches 1. By symmetry, the curvature also accumulates at the 6 other shrinking 
rational curves, {z l = C} and {z 2 = C} for C = |, |, \ + |. 

In the z % coordinate system used in figures || and ^, the exceptional divisor (of the original 
orbifold) is represented as a single point, namely the origin. This is misleading, since in fact 
it is topologically an S 2 . It's interesting to study how its geometry changes as we vary a. 
The induced metric on the exceptional divisor of the Eguchi-Hanson geometry is that of a 
round sphere. We therefore expect the same to be true for the exceptional divisor of the 
Kummer surface at small a. In figure [3|, we plot the Ricci scalar of the induced metric on the 
exceptional divisor. For small values of a, it is indeed uniform. However, as a increases, it 
becomes non-uniform (its integral is of course fixed at 8ir): the sphere is becoming prolate. 
The poles, where the curvature is highest, are at y = and y' = 0, in other words where the 
exceptional divisor intersects the planes z 1 = and z 2 = respectively; these are the points 
that are closest to neighboring exceptional divisors. 

Similarly, it is interesting to study the induced geometry of the shrinking rational curves, 
for example the curve {z 1 = 0}, as shown in figure ||. At a = its geometry is that of a flat 

4 Animations of these plots as a function of Im z 2 , showing the entire four-dimensional geometry, are available 



for download at http://schwinger.harvard.edu/~wiseman/k3/ 
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Figure 1: Isosurfaces of the Euler density p in the z l coordinate system, on the Imz 2 = slice. 
Only the fundamental domain of the discrete symmetry group is shown, namely the cubical re- 
gion < Re z 1 , Imz 1 , Re z 2 < \ (thus the full Imz 2 = slice, which is a T 3 /Z 2 , is 32 times 
as large as what is shown). Blue indicates lower curvature, red higher, with surfaces shown at 
p = 1CT 1 - 5 , KT 1 , . . . , 10 4 , 10 4 - 5 . Left to right and top to bottom, a = 0.03, 0.13, 0.28, 0.50. 

T 2 orbifolded by Z2, which is topologically S 2 but with all the curvature concentrated at the 
4 fixed points z 2 = 0, ^, |, \ + \ (picture a square envelope). As a increases, the curvature 
spreads out around the sphere, which becomes rounder and rounder. In the limit a. — ► 1, the 
geometry in the vicinity of the shrinking rational curve should approach that of an Eguchi- 
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Figure 2: As in figure but for a = 0.61, 0.72, 0.85, 0.92. 



Hanson metric whose exceptional divisor has area A = \b 2 - 2vra 2 . Indeed, we see that for 
values of a approaching 1, the sphere becomes almost completely round. As another test that 
the geometry is approaching Eguchi-Hanson, in figure 4 we plot the maximum and minimum 
values of p on the curve against a. On an Eguchi-Hanson, the Euler density is constant on 
the exceptional divisor; the solid curve in that plot is the value of p on the exceptional divisor 
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12 3 



Figure 3: Log (base 10) of the Ricci scalar of the induced metric on the exceptional divisor w = 0. 
Points on the sphere are mapped stereographically from their y coordinate, with the north and south 
poles corresponding to y — and y' — 0. Values of a are as in figures |l| and |^. The average value of 
the Ricci scalar over the sphere decreases as its area increases, as is clear from the figure. 

of an Eguchi-Hanson of the appropriate size. One sees that as a approaches 1 the maximum 
and minimum values of p approach each other and the Eguchi-Hanson value. 

Given the Ricci-flat metric, one can compute the spectrum of various geometric operators 
of physical interest, the simplest example being the scalar Laplacian. As a proof of principle, 
we used our numerical metrics to compute a low-lying eigenvalue (and eigenfunction) of it. 
The eigenfunctions of the Laplacian on our symmetric Kummer surface can be classified by 
their eigenvalues under the Z| translation subgroup of its full symmetry group. We calculated 
the lowest eigenvalue in the sector with eigenvalue —1 under all 4 translations. The results 
are shown in figure 5. The eigenvalue at the orbifold point a = 0, which can be computed 
exactly, is also shown for comparison. The fact that the eigenvalue increases with a can be 
understood heuristically (at least for small values of a) in the following way. As the rational 
curves discussed above shrink, the geodesic distance between a point on the exceptional 
divisor, such as y = w = 0, and its image under each of the discrete translations decreases, 
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12 3 



Figure 4: Log (base 10) of the Ricci scalar of the induced metric on the rational curve {z 1 = 0}. (The 
Ricci scalar is everywhere positive.) Using the function y(z 2 ) = —ip(z 2 )/p(^) (with p the appropriate 
Weierstrass elliptic function), points are mapped from T 2 /Z2 to the complex plane; y is then mapped 
stereographically to the sphere. The maxima of the Ricci scalar are at z 2 = 0, i, |, h + 5. Values of 
a are as in figures g a nd |. The average value of the Ricci scalar over the sphere increases as its area 
shrinks (see equation ( |3.16| )), as is clear from the figure. 

requiring steeper gradients in the eigenfunction. Finally, we would like to point out that, as 
we are studying a long wavelength eigenfunction, even the lowest resolution (A) computes 
the eigenvalues to typically within a few percent of the continuum extrapolated value. 

We close this section on a technical note. Careful inspection of the isosurfaces plotted for 
the smallest value of a in figure [l] shows a localized deviation from sphericity, caused by small 
errors near patch overlaps. Similarly, in figure |3| the smallest a sphere is slightly less uniform 
than for the next larger value of a; in figure || for the largest value of a we can just see some 
"rings" where small errors are introduced due to coordinate patch overlaps; and the distance 
from the extrapolated continuum values in figure |6| increases near a = or 1. Numerical 
discretization errors are to be expected, and will be larger where there are higher curvatures. 
Our solutions containing the highest curvature regions occur near the two orbifold points 
a = 0, 1, explaining why we see the effects mentioned above. Hence, at a fixed resolution, 
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Figure 5: Maximum and minimum values of the Euler density p on the rational curve {z 1 = 0} as 
a function of a. The solid curve is the value of p on the exceptional divisor of an Eguchi-Hanson 
geometry whose exceptional divisor has area A. 



the global quality of the solutions will not be as high near these orbifold points, although 
the local geometry away from the regions of high curvature should be quite acceptable. Of 
course, increasing the resolution, the quality of the solution will improve, independently of 



where we are in moduli space. These issues are discussed in more detail in appendix A. 3 



4. Discussion 



4.1 Generalizations 



We opened this paper with the question of whether it is possible, in practice, to solve the 
Einstein equation numerically on a Calabi-Yau manifold. We have shown that the answer is 
affirmative when the Calabi-Yau is a K3 surface with a high degree of discrete symmetry. In 
this subsection we will investigate the possibility of generalizing this accomplishment, first 
to more general K3 surfaces, then to Calabi-Yau three-folds, and finally to Kahler manifolds 
with cosmological constant or with matter. 

It is clear that in principle our method extends to a general blow-up of a Kummer surface, 
and more generally to any K3 surface. It is merely necessary to implement the topology and 
complex structure with complex coordinates defined appropriately on an atlas of charts — for 
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Figure 6: An eigenvalue of minus the scalar Laplacian as a function of a. Specifically, the lowest 
eigenvalue whose eigenfunction is antisymmetric under all 4 generators of the Z| translation group of 
the symmetric Kummer surface. A, B, and C represent increasing lattice resolutions, as detailed in 
Appendix The exact eigenvalue at the orbifold point a — is also shown. 

example one could use patches derived from any algebraic construction of K3. However, 
the question remains as to what resolution is attainable in the absence of large amounts of 
discrete symmetry. Our highest resolution, D, simulated approximately 2 x 10 7 ~ 80 4 points. 
However, due to the high degree of discrete symmetry our one parameter family of K3's enjoy, 
every computed point actually represents 2 9 points in the true K3 geometry, and hence we 
have described the full K3 with an effective resolution of around 400 4 . 

On a current high-end desktop computer (with 1 to 2 gigabytes of memory) one could 
comfortably increase the total number of points simulated to 10 8 . For a K3 with no discrete 
symmetry this would yield a resolution of 100 for the full K3 geometry. To estimate how 
accuately this could represent the metric, we may compare it to our intermediate resolution, 
B. With a linear resolution 4 times lower than that for D, the effective resolution for the full 
geometry is also about 100 4 . As seen in subsection |3.2j and appendix^ we find that resolution 
B adequately reproduces the geometry and derived properties, such as the low wavelength 
eigenmodes of the Laplace operator, provided one is not too close to the edge of the Kahler 
cone. For example, in figure |6| one sees that run B computed the Laplacian eigenvalue with an 
accuracy of around 1% compared to the extrapolated continuum value. Near the edge of the 
Kahler cone, where regions of high curvature develop in the manifold, the best strategy may 
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be to combine numerical with analytic techniques. For example, at values of the moduli near 
an orbifold point, one could patch an analytic Ricci-flat metric, such as the Eguchi-Hanson 
metric, into the region where the high curvature is developing. 

If we now wish to move to a Calabi-Yau three-fold, 10 s points translates to a mere 20 
points linearly in each direction. This is over a factor of 2 less in linear resolution than 
the lowest effective resolution used in this work (run A, which had an effective resolution of 
50 4 for the full geometry). This might be acceptable if one were well away from the Kahler 
cone edge, but is certainly rather low. On the other hand, if one were to consider a highly 
symmetric three-fold, such as a symmetric blow-up of T 6 /Z3, then one would again expect 
to attain an effective resolution of around 100 points linearly in the full geometry, and thus 
again expect accuracy comparable to, or better than, our resolution B. 

Moving to the general three-fold appears to be a very challenging task. As we discuss 
in Appendix [Aj processing time actually scales rather well with increasing dimension. In- 
stead the problem is limited by storage. In six real dimensions any appreciable increase in 
linear resolution is extremely costly. Thus whilst 20 6 points would be possible on a desktop 
computer, 40 6 requires 64 times more memory and is already beyond the abilities of modest 
clusters. Often a tough computational problem becomes easy in time, as computer memory 
and speed have closely followed Moore's prediction of a doubling every 2 years. However, even 
assuming that Moore's law continues to hold, it will require 12 years to increase the linear 
resolution by a factor of 2 for the three-folds. It is therefore clear that in order to tackle the 
general three-fold, one must employ considerably more sophisticated discretization schemes 
than we have used here, presumably adapting the lattice points to regions of high curvature. 
Whilst adaptive grids are difficult to implement in the elliptic context, one can easily use 
fixed grids that increase resolution in areas where curvature is expected to be high. 

To summarize, we expect that for a general K3 surface one can obtain very satisfactory 
results provided one does not wish to probe too near the edge of the Kahler cone. For a 
highly symmetric three-fold similarly high quality results can be expected. However, moving 
to the general three-fold appears tough and new techniques must certainly be employed. 

Having considered the problem of constructing K3's at arbitrary points in moduli space, 
we should point out the enormity of the moduli space itself: with 57 directions to explore it 
is highly implausible that one could ever map the entire space of K3's. On the other hand 
it is unlikely that one would ever need to map the entire space. One can imagine wishing to 
find K3 surfaces with specific properties; presuming the observables of interest vary smoothly 
over the moduli space, it is plausible that one could scan the moduli space for examples that 
fit the specific requirements. 

Whilst Calabi-Yau's have large numbers of moduli, spaces with Ricci curvature tend to 
have fewer of them. An example of this are the four (real) dimensional del Pezzo surfaces 
dP n (n = 1, . . . , 8), compact Kahler spaces that can be constructed by blowing up n points in 
CP 2 . Del Pezzos have positive first Chern class, and those with n > 3 admit Kahler-Einstein 
metrics with no Kahler moduli. The cases of dP3 and dP4 are particularly interesting because 
they also have no complex structure moduli. The Einstein equation can again be reduced 
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to a Monge-Ampere equation for the Kahler potential. As in our K3 example, two patches 
would be required to cover each of the n blown up points, and three to cover the ambient 
CP 2 . Hence we expect the methods we have applied here for K3 to be applicable, and the 
Kahler-Einstein metrics to be attainable to high resolution on a desktop computer (certainly 
comparable to or better than our resolution C, depending on how many points are blown up 
in CP 2 ). This would be very satisfying for dP3 and dP4 as, with no moduli, one would in 
principle have constructed these geometries explicitly and completely. 

Other physically interesting geometries that are related to Calabi-Yau manifolds are 
supersymmetric flux compactifications in string theory. A class of type IIB solutions can be 
constructed as warped products of flat four-dimensional Minkowski spacetime and a Ricci-flat 
Calabi-Yau three-fold, where the warp factor satisfies a Poisson equation on the Calabi-Yau 
sourced by fluxes, D-branes, and orientifold planes (see e.g. [19, 20, ^2|). As we discussed 
above, finding the metric on a generic three-fold is probably out of reach using the methods of 
this paper. However, considering fluxes on K3 or a highly symmetric three-fold would likely 
be a manageable task. Having found the Ricci-flat metric, solving the Poisson equation on 
this geometry is quite simple — indeed even easier than finding eigenfunctions of the Laplacian 
as we did earlier in this paper. Studying the solutions to the Poisson equation on the Calabi- 
Yau background would provide a detailed understanding of how the fluxes backreact on the 
vacuum geometry, and in particular of how fluxes on adjacent cycles interact. 



4.2 Lessons for solving general Euclidean geometries 

The key simplification in our work has been that of Kahler geometry. What are the prospects 
for constructing general Euclidean geometries? Without Kahlerity we require many metric 
components to describe the geometry, and the Einstein equation becomes complicated. Whilst 
this may be technically complicated, in principal one would hope that using a harmonic gauge 
condition would allow the system to be locally solved as an elliptic relaxation problem. Note 
that one could also use a gauge fixed Ricci flow, but as discussed eariler, it is more efficient 
to solve the elliptic Ricci flatness condition directly rather than to construct an entire flow 
when only the endpoint is required. In our K3 example, the complex coordinates on the 
coordinate patches provide exactly such a local harmonic set of coordinates. However the 
most challenging aspect of the problem is to understand the global issues, such as residual 
coordinate freedom, adaptedness of the coordinates, and moduli of the solutions. 

Let us for a moment consider the problem of finding the Ricci flat metric of symmetric 
K3's as we have done, but ignoring the Kahler structure and using only real geometry. We 
might hope to find harmonic coordinates on the various patches that make up the topol- 
ogy. This would require us to solve the harmonic gauge condition (essentially locally solving 
Laplace equations) at the same time as the Einstein equation. Presumably this full system 
is globally elliptic in the case of K3, although for more general geometries we should note 
that negative modes of the Lichnerowicz operator will exist, as occurs for the Euclidean 
Schwarzschild solution, and it is unclear how these would affect the situation. 
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Given the link between the complex coordinates of Kahler geometry and the harmonic 
coordinates natural for Euclidean real geometry, we can make various speculations. We saw 
that our complex coordinates were well adapted to the symmetric K3 geometry, and one 
might hope the same to be true for more general harmonic coordinates. As explained above, 
for the complex coordinates there are no residual holomorphic coordinate transformations; 
in real geometry, choosing harmonic coordinates one again expects only finitely many resid- 
ual coordinate freedoms on a compact manifold. Assuming this, one might conclude from 
our work that we may see the complex structure moduli of the K3 arise simply from the 
global data required to specify the harmonic coordinates, just as we have fixed the complex 
structure moduli by taking particular complex coordinates on the manifold. Then for more 
general compact manifolds, one might plausibly associate physical moduli to global choices 
when constructing the harmonic coordinates (that are not simply one of the finite residual 
coordinate transformations). 

Clearly about any Ricci flat real geometry one can always linearize metric fluctuations and 
then, in principle, directly determine the zero modes of the resulting Lichnerowicz operator, 
and hence determine all physical moduli of the solution. However, this is obviously very 
complicated to imagine doing in practice, and what we really wish to find is a way to include 
the moduli as boundary conditions in the elliptic problem as we have done in our Kahler 
example. Assuming our presumptions above about the complex structure moduli hold, the 
key remaining question is how to understand the Kahler moduli of K3 as boundary conditions, 
without actually making use of the Kahler structure. It might be possible to understand this 
in terms of the volumes of minimal representatives of the two-cycles. However, whilst for K3 
this approach would work, for the Kahler-Einstein case of the del Pezzos it cannot, since the 
two-cycles present are not associated with any moduli; hence this approach would not work 
in general. Thus, finding new ways to understand the Kahler moduli and how to actually 
implement fixing them whilst solving the real geometry Einstein equations for K3, look to be 
important questions. If they can be addressed, it might allow one to understand the moduli 
of general real geometries, and enable explicit metrics to be found in very general Euclidean 
geometry-matter systems. 

4.3 Applications 

We now briefly discuss possible applications for the numerical construction of geometries. We 
consider first mathematical, then formal physical, and finally phenomenological applications. 

There are various outstanding mathematical conjectures regarding the geometry of Calabi- 
Yau manifolds. Most notorious is the Strominger-Yau-Zaslow conjecture, that Calabi-Yaus 
with mirrors can be constructed as toric fibrations and that mirror symmetry acts by T- 
duality on the fibers 24|. A related conjecture concerns mean curvature flows and special 



Lagrangian submanifolds in Calabi-Yau manifolds (25|]. In principle these conjectures can be 
tested directly on any Calabi-Yau for which the Ricci-flat metric is known explicitly. 

From the point of view of physics, explicit constructions of Ricci-flat Calabi-Yau metrics 
may help us learn more about the sigma model description of geometry. The Ricci-flat metric 
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on K3 is the target space of an Af = (4, 4) non-linear sigma model. Due to the high degree 
of supersymmetry the classical Ricci flat metric receives no perturbative or non-perturbative 
corrections in ol fl2|| . Thus the geometries we have constructed in this paper can, remarkably, 
be viewed as fully quantum geometries from the viewpoint of this sigma model. Knowing these 
metrics then in principle allows one to compute properties of the quantum sigma model — for 
example, the spectra of operators on the target manifold correspond to conformal weights on 
the worldsheet. 

The Ricci-flat Calabi-Yau three-folds are again target spaces of sigma models. However, 
now the classical geometry only gives the leading a', or large volume, approximation to the 
true quantum geometry. Understanding how a' corrections modify the classical geometry 
is an important physical issue. Supersymmetry implies that these corrections preserve the 
Kahlerity of the metric, and hence they will appear as higher derivative terms modifying the 
Monge- Ampere equation for the Kahler potential. Each higher derivative term will make this 
equation less local, but in principle we may still apply the same local iterative methods we 
have used here to solve it. (Presumably on a compact manifold, one could in principle include 
infinitely high derivative terms if their form were known.) From a discretized viewpoint, if we 
linearized the equation about some background, we would find anJVxJV operator (N being 
the total number of lattice points) and the structure of its matrix representation will no longer 
be sparse. However, the terms that fill in the zero components in the sparse Monge- Ampere 
case will be small, being down by factors of a' . Hence our Gauss-Seidel iterative methods 
may still work, although obviously evaluation of the equation at each point will take much 
longer. 

Finally, from a phenomenological point of view, being able to compute metrics is crucial 
for actually making contact with low energy physics in string theory. Whilst for simple 
vacuum Calabi-Yau reductions it is possible to compute the entire low energy effective action 
using only topological data [26|, as soon as matter is added to the compactification manifold 



this is no longer true. The simplest example is adding a single brane that is localized in the 
compact space. The moduli space metric for this brane, and hence the kinetic term for its 



position in the low-energy action, is simply given by the metric on the internal space [27]. 
Thus from a physics standpoint, being able to compute the low energy action of geometric 
reductions is a strong motivation to further understand and improve the numerical geometry 
methods we are exploring here. 

It is worth mentioning that in TeV fundamental scale senarios |28| . |2?| it is possible that 
in a few years the LHC might directly probe not just the low energy action, but also high 
energy excitations on the internal space, such as Kaluza-Klein modes. Although the LHC 
would likely not measure a sufficient number for any detailed spectroscopy of the internal 
space, future colliders would then be able to measure these massive excitations accurately. 
If one could extend our methods to the general three-fold — and surely this would be enough 
motivation to direct serious resources to it — then one might scan the moduli space of Calabi- 
Yaus, presumably with fluxes, to find candidate geometries matching the observed resonances. 
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A. Details of numerical construction 

A.l Construction of the atlas and initial data 



As discussed in section |3.l| our one parameter family of K3's have many discrete symmetries. 
In particular these imply we may take the Kahler potential to be identical in each of the 16 
regions describing the blow up of the torus fixed points. These regions are each described 
in terms of 2 coordinate patches given by w,y and w',y' and the symmetry z 1 <-> z 2 implies 
that the Kahler potential may be taken to be identical in each of these. Thus we reduce our 
problem to one patch describing the fundamental domain of the torus using z 1 , z 2 coordinates, 
where we orbifold by the identification (z l ,z 2 ) = (— z 1 ,— z 2 ), and one patch describing (half 
of) the blow up of the fixed point contained in that fundamental domain using coordinates 
w, y. We term these patches the 'Torus patch' and the 'Eguchi-Hanson patch'. To begin 
describing the geometry of our atlas we define, 



a = \z 1 \ 2 + \z 2 \ 2 



2M(i + M 2 ) (A.i) 



We take our torus patch to cover the coordinate range of the fundamental domain, but exclude 
the region near the blown up fixed point so, 

1 11 1 

— <Rez ll2 <- , — <Imz 1,2 <- and a > a min (A.2) 
4 4 4 4 

Outside the boundaries of this domain we must act with the translations z 1,2 — ► z 1,2 ± (1, i)/2 
to map the point back into the domain. Note that when we do this, we must also ensure 



we perform the torus Kahler transformation derived from 3.11 reduced to this domain. The 
Eguchi-Hanson patch is taken to have coordinate range, 

< \y\ < y m ax and a < a max (A.3) 

In order to cover the manifold we must ensure y m ax > 1, and a max > a m i n with a max < 1/4? 
to avoid complicated multiple overlaps. 

In order to evaluate our Monge- Ampere equation at the edge of a coordinate patch we 
must necessarily compute derivatives involving Kahler potential data from neighbouring co- 
ordinate patches. Once we have finite differenced our patches this will require the Kahler 
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potential at some point to be computed from the neighbouring patch and then Kahler trans- 
formed into the original patch. The coordinate location in the neighbouring patch will not 
necessarily fall on a lattice point and therefore we will need to perform interpolation to com- 
pute the desired Kahler potential 5 . Thus we require our patches to overlap sufficiently in 
order to perform our necessary interpolations. 

For example, suppose when we evaluate a derivative at the boundary of the Eguchi- 
Hanson patch we require knowing the Kahler potential at a point still with a < a max but 
now \y\ > y m ax- Then we must use the coordinate patch w',y' - which due to the discrete 
symmetries is identical to the w, y one. Explicitly, we transform to the w', y' coordinates where 
still <r < a m ax, but now \y'\ < y m ax, so the point does indeed lie within this w',y' patch. We 
find the Kahler potential in this patch using the necessary interpolation, and then return to 



our original coordinate patch w, y by performing the necessary Kahler transformation |3.14 . 
Similarly, at the large a boundary of the Eguchi-Hanson patch, or small a boundary of the 
torus patch we will find the Kahler potential in that patch by interpolating from the other 
and performing the appropriate Kahler transformations 3.12| , |3.13| . 



We note that in fact we actually require only one coordinate patch, as the Eguchi-Hanson 
patch can quite satisfactorily represent the torus region. However, the torus patch boundary 
conditions, essentially derived from |378 - p . 11 become complicated and non-local in the w,y 



coordinates, and we have found it simpler to use the two patches above, rather than the 
minimal choice of one patch. 

Even after the reduction to these 2 patches, and the coordinate domains above, we still 
have discrete symmetries remaining. The orbifold symmetry, holomophic isometries z 3 — * iz 3 , 



z 



1 z 2 and anti-holomorphic isometry z ' 2 <-> z 1,2 further reduce the torus coordinate 



domain (in addition to a > a m in) to, 



0<Rez 1,2 <- , 0<lmz 1 ' 2 <- 
- 4 - 4 



and (z\z 2 ) ~ (z 2 ,z r ) 

(z\z 2 ) ~ (z\z 2 ) (A.4) 



In the Eguchi-Hanson patch these isometries allow us to further reduce A. 3 to, 



0<Rew,y <Imw,y (A. 5) 

The Kahler potential does not transform under these discrete isometries, and therefore if we 
now require a point to be evaluated outside our reduced coordinate domains we simply use 
these discrete symmetries to map the point back into our reduced coordinate domain. 

In the data we present we have chosen a fixed geometry for the patches and their overlaps. 
This is independent of the numerical resolution of the discretization so that for convergence 

5 Note that whilst extrapolation is more economical as we would require no patch overlap, it is also po- 
tentially dangerous as it is likely to introduce spurious data, and since specifiying the correct data is our key 
concern we have opted to use interpolation which takes a little more storage (due to the patches overlapping) 
but removes the risk of specifying data incorrectly. 
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testing we are comparing like with like. As discussed in section 2.3 it also allows us to directly 
compare the Kahler potential in the same overlapping regions as we vary the resolution which 
provides a check of numerical convergence. The parameters we have chosen are, 

Vmax = 1-25 a min = -j x 0.32 and a max = X 0.60 (A. 6) 

Before we begin relaxing the Monge- Ampere equation we require some smooth initial data 
compatible with our choice of Kahler parameters a, b (although in fact we have found that in 
some cases even taking non-smooth initial data the algorithm still converges). We construct 
this by taking an initial guess Kahler potential to have the behaviour of the Eguchi-Hanson 
potential near a = in the Eguchi-Hanson patch, and which then interpolates smoothly up 
to second derivatives to behave as the flat torus potential for a > a max . 

The Eguchi-Hanson potential in the coordinate patch w, y is given by, 

1 + v 1 + ^ 

and the flat torus potential is simply K tOTUS = \ob 2 . For our interpolation, in our Eguchi- 
Hansen patch (with a < o max ) we take, 

K (y,w) = ^a 2 log (1 + \y\ 2 ) + k + k 2 a 2 + k 4 a 4 

with k = — (6a 2 + 3a max b 2 - 8a 2 log a max ) 
lb 

k2 = — (3a m axb 2 - 4a 2 ) 

""max 

k4 = -TjT-g (Vmaxb 2 ~ 2a 2 ) 

(A.8) 

and in the torus patch with a > a m i n we take, 

G < Vmax K(zl,z2) = ^ lo g a + k o + k 2 (J 2 + A; 4 cj 4 

1 2 

O" > Vmax K(z\,z2) = 2 

(A.9) 



' max • 



The constants above ensure the Kahler potential has smooth second derivatives at a = ay, 
However, it is easy to show that the above does not define a positive Kahler form over the 
whole range a = to 1. As discussed in the main text, one advantage in solving the Monge- 
Ampere equation directly, rather than performing Ricci flow, is that the Kahler form need 
not be positive. Hence the simple interpolation above suffices as initial data. 
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A. 2 Discretization, memory and time requirements 

We discretize our system in the most naive way. In both patches we discretize by creating 
a uniform lattice in the real and imaginary parts of each complex coordinate. We then use 
second order finite differencing to implement the Monge-Ampere equation and third order 
accurate interpolation at the patch overlaps. 

At each step of our relaxation, we firstly interpolate the values of the Kahler potential at 
the very edges of a patch from the appropriate neighbouring patches, and secondly perform one 
iteration of the Gauss-Seidel generalizated to our Monge-Ampere equation. When updating 
each point, as a by-product we may quickly compute the value of det g at that point. During 
the Gauss-Siedel update, we keep a running average of this quantity over all lattice points, 
and then use this as the value of A for the next step. This ensures our Gauss-Seidel iterations 
asymptote to a fixed solution. 

In this paper we have used 4 different resolutions each differing in linear resolution by a 
factor of 2. In the the torus patch we discretize with equal lattice spacing dz in the real and 
imaginary z l and z 2 directions. In the Eguchi-Hanson patch we discretize with spacing dw 
in the real/imaginary w directions and dy in the real/imaginary y directions. We label the 
various resolutions A — D, and they are specified as, 



(A.10) 



Remembering that the coordinates z range from — > 0.25 for our reduced domain in the 
torus patch, this yields 80 points along a side of the torus in resolution D. 



A 


dz = 


0.025 


dw = 


0.003125 


dy = 


0.25 


B 


dz = 


±0.025 


dw = 


±0.003125 


dy = 


±0.25 


C 


dz = 


-^0.025 


dw = 


^0.003125 


dy = 


^0.25 


D 


dz = 


^0.025 


dw = 


2^0.003125 


dy = 


^0.25 



In the torus patch we implement the 2 identifications in A. 4 imperfectly by simply only 
storing points with Im(z 2 ) > max (Re(z 1 ), Irn^ 1 )) . This is rather convenient, but doesn't 
fully take advantage of these discrete symmetries. Asymptotically this yields a reduction of 3 
rather than the optimal value 4 for these 2 identifications. With this minor imperfections in 
mind, the total number of points stored in computer memory for each resolution to represent 
the manifold is then, 

number of points A = 6 x 10 3 

B = 7 x 104 , a ■, ■, s 

« (A.ll) 
C = 1 x 10 6 v ' 

D = 2 x 10 7 

The Gauss-Seidel scheme is implemented by solving the discrete Monge-Ampere equations 
at each point in the lattice. We found that under-relaxation was not required for stability 
and the scheme converged stably. Interestingly we found that we could not over-relax the 
equation with any appreciable over-relaxation rendering the scheme unstable. Presumably 
this is an effect associated with the patch boundaries rather than their interiors which we 
expect to behave in an analogous manner to the Poisson equation. Thus in principle it might 



- 28 - 



be possible to include an over-relaxation parameter that varied over the patch to be one at 
the edges, but larger than one in the interior. 

We measure the distance from convergence by computing the maximum update of the 
Kahler potential in a given Gauss-Seidel iteration over the whole lattice. This is equivalent 
to computing the maximum violation of the discretized Monge- Ampere equation. 

When this number falls below 10 -12 we classify the solution as having relaxed. Certainly 
for any quantity we have computed, there is no further change if the solutions are subjected 
to further iterations of the Gauss-Seidel scheme. The time taken for our implementation of 
the algorithm from initial guess to the relaxed condition is shown in figure [7] for the various 
resolutions A — D on a standard desktop computer (3Ghz Pentium, 500 Mb). After a little 
relaxation, as for the Poisson equation, the number of Gauss-Seidel iterations required to 
improve the discretized error by a factor of 10 quickly tends to a constant. This is also shown 
in the same figure. 
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Figure 7: The lefthand plot shows the log (base 10) total time measured in days for the various 
resolutions A — D to converge to our required degree. The righthand plot shows the number of Gauss- 
Seidel iterations to improve the maximum error in the discretized Monge- Ampere equation by a factor 
of 10. 

We see that both these quantities exhibit only a weak dependence on the position in 
moduli space that we choose. The lowest resolution A relaxes in seconds while our highest 
resolution D requires up to a week. 

The time taken to relax using the local Gauss-Seidel scheme can be estimated as going 
as N 1+2 / d in d real dimensions where iV is the total number of lattice points. Every iteration 
takes a time of order N, and the total number of steps can be estimated as N 2 / d by consid- 
ering the spectral radius of the linearized Laplace operator. We see this scaling is certainly 
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consistent with the results of figure ^. 

At every step in the iteration we must also perform an interpolation of points at the edge of 
each coordinate patch from its neighbours. We only require the edge points to be interpolated 
that are required by our second order differenced Monge-Ampere equation. Thus the work 
involved scales as a codimension one quantity, namely as In practice our third order 

interpolation actually takes considerable time. Asymptotically at large N, however, it will 
obviously become subdominant to the Gauss-Seidel iteration time which scales as N. 

It is very interesting to note that as we move to higher dimensions, the total relaxation 
time more and more closely approaches N. Thus the advantage in using highly non-local 
schemes such as multi-grid to improve convergence times (typically to iVlogiV) becomes 
considerably reduced. In this sense we may claim that the problem of extending our methods 
to Calabi-Yau 3-folds is storage limited rather than speed limited. 

In a memory limited problem, it is rather natural to move to higher order methods. 
Our implementation uses second order finite differencing on the Monge-Ampere equation. 
However, we might hope for improved convergence to the continuum if we were to use 4 th order 
discretization of the Monge-Ampere equation (and also for interpolating between patches). 
We did try this in our case of K3, but the additional time required by each iteration slowed the 
total convergence time to approximately the same as the next higher resolution using second 
order differencing. In order to precede to our highest resolution we decided to stay with 
second order differencing. Bear in mind that whilst a higher order method approaches the 
continuum more quickly, in order to resolve short length scales one requires high resolutions, 
so to get accurate results near the orbifold regions of our moduli space, we require the highest 
resolutions possible. 

With a particular physics or maths question in mind, one might attempt computations 
using more resources and improved stamina than we have used here. Then the fundamental 
problem of limited storage should probably be tackled by both a combination of more efficient 
discretization, and also higher order methods. The increased cost in processor time could 
certainly be ameliorated by some form of parallelisation - which is well suited to this problem 
which, afterall, naturally divides into coordinate patches. 

The eigenfunction of the scalar Laplace operator presented in section |3.2| was computed 
using a naive iterative scheme. The initial guess for the eigenfunction if)( a \ with odd parity 
under the Z\ translation isometry {z l —> z l ± 1/2) was taken to be that for the torus orbifold, 
a = 0, where, 

^(o) = cos 2ir Re z 1 cos 2tt Im z 1 cos 2tt Re z 2 cos 2tt Im z 2 (A. 12) 

giving an eigenvalue 4(27r/6) 4 . The eigenvalue equation was solved using Gauss-Seidel itera- 
tion everywhere but at one point, with the eigenvalue being determined dynamically from the 
condition that the eigenfunction be smooth at that point. Of course, explicit independence of 
the actual point chosen was checked. This method is very simple to implement, but is really 
only suited to finding the lowest eigenfunction in a parity sector. More sophisticated (but 
standard) methods would be required to computer higher eigenfunctions. 
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A. 3 Convergence tests 

We now breifly present data that demonstrates increasing numerical resolution improves ob- 
served quantities in a manner consistent with a second order approach to the continuum. 

Firstly as discussed in the main text, we compute a numerical A rather than using the 
analytic value A. This ensures the numerical Monge- Ampere equations converges to a static 
end point. We may then compare this numerical determination of A with the true analytic 
value. In figure ^ we plot the difference between these values for the various resolutions as 
a function of position in moduli space. We see the error is greatest near the orbifold point 
at a. = 1 as expected. We clearly see that the errors improve with increasing resolution 
consistent with second order scaling - we remind the reader that the 4 resolutions A — D each 
differ by a linear resolution factor of 2. 
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Figure 8: Plot of the log (base 10) absolute difference between the numerically determined A and the 
true analytic value A for the various resolutions A — D. 

In figure || we plot the integrated Euler density for the resolutions B — D. The result, the 
Euler number, has true value 24. We clearly see here that increasing resolution does indeed 
improve the numerical determination of this quantity as we would hope for. The lowest 
resolution gives a very poor estimate of the Euler number and we have not included it here. 
Resolution B still provides a rather poor estimate. We see the error in the Euler number is 
practically quite small for resolution C provided we are not too near either orbifold point, 
and resolution D gives errors of less than 0.1% over most of the moduli space. 

Finally our patches geometrically overlap in coordinate space. Therefore as discussed 
in the main text, once we have found a solution, we can test its quality by computing the 
error in the Kahler potential in the overlapping regions (obviously taking into account the 
relevent Kahler transformation between the 2 patches). In figure [k] we show the maximum 
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Figure 9: The Euler number integrated from the Euler density for the resolutions B, C, D. 



error found by comparing 2 different patch overlaps; the overlap of the Eguchi-Hanson patch 
with itself, and then its overlap with the torus patch. We see that this maximum error again 
decreases consistent with second order scaling as resolution is increased. 



B. Homology of Kummer surfaces 

In this appendix we derive some properties of the second homology group of Kummer surfaces 
that were used in the main text. We also show how to relate the natural basis for the homology 
to the standard basis for the integral homology of K3. As far as we know this explicit relation 
is new. 
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Figure 10: Lcfthand plot shows the log (base 10) absolute maximum error between Kahler potentials 
on the overlap of the Eguchi-Hanson patch with itself. The righthand plot shows the maximum error 
on the overlap between the Eguchi-Hanson and torus patches. 



K3 has second Betti number 62 = 22. In the Kummer construction, 6 two-cycles are 
inherited from T 4 , and the other 16 are the exceptional divisors at the blown up fixed points. 
The second homology of T 4 is generated by the 2 holomorphic curves {z 1 = C} and {z 2 = C}, 
as well as the 4 non-holomorphic curves {Re z 1 = C , Re z 2 = C 2 }, {Re z 1 = C , Im z 2 = C 2 }, 
{Imz 1 = C^Rez 2 = C 2 }, and {Imz 1 = C^Imz 2 = C 2 }. When taking the orbifold, one 
may choose the constant (s) in a such a way that the curve either passes through or avoids the 
fixed points (e.g. the holomorphic curves pass through four fixed points if C = but avoids 
them if C = 1/4). In the latter case, one must include the image under the orbifold, e.g. 
{z 1 = C} U {z 1 = —C}, to obtain a two-cycle in the orbifold. (The cycles that pass through 
the fixed points are linear combinations of those that don't and the exceptional divisors.) For 
our purposes, it is simplest to take as a basis the 6 cycles that avoid the fixed points, which 
we will refer to as "torus cycles", along with the 16 exceptional divisors. We refer to these 
cycles as ci, where c± 2 are the holomorphic torus cycles, 03^5 6 are the non-holomorphic torus 
cycles, and 07^.^2 are the exceptional divisors (which are holomorphic). 

It is straightforward to write down the intersection matrix in this basis. The torus cycles 
intersect each other in pairs, with intersection number 2. For example c\ intersects C2 at two 
points, (C, C), and (C, — C) (of course we don't count (— C, — C) and (— C, C) separately). 
By construction, none of the torus cycles intersect the exceptional divisors. The latter do not 
intersect each other, but have self-intersection —2 (since they are topologically Ci-^'s). All 
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in all, we find the following block diagonal intersection matrix: 

~U 



hi j = #(c/,cj) = 2 



U 



u 



'16 



u 



1 

1 



(B.l) 



From the Kahler transformations (3.8- 3.14 ) the periods of the Kahler form jj = J ci J 
may be computed for the symmetric Kummer surfaces: 

V, 7 = 1,2 

31 = \ 0, 1 = 3,4,5,6 . (B.2) 
na 2 , 1 = 7, ...,22 

In the case of the holomorphic curves (I = 1,2,7, .. . ,22) these periods are their areas. In 
terms of the periods the volume is 



1 



1 



V = ~ / J A J = ^h IJ j!jj = ^ - 4^a 4 , 



1. 



,2 4 



(B.3) 



where h IJ is the inverse intersection matrix, 



JJ 



U 



U 



u 



-I 



16 



(B.4) 



Consider now the holomorphic curve c\ = {z 1 = 0}. This intersects C2 at one point, 
(0,C), but none of the other torus cycles. It also intersects, at one point each, 4 of the 
exceptional divisors, namely those located at (0, 0), (0, ^), (0, |), (0, \ + |), which we will call 
c 7,8,9,io- Using the intersection matrix, we have 



ci = - (ci - c 7 - c 8 - c 9 - Clo) 



Hence the area of this curve is 



A 



Cl 



2ira z , 



(B.5) 
(B.6) 



as claimed in subsection 3.1. 

Now let us return to the general Kummer surface. The two-cycles cj are obviously 
integral. However, since the inverse intersection matrix ( |B.4| ) is not composed of integers, 
they do not form a basis for the integral homology H2 (K3, Z), but only a sublattice of it. The 
standard basis {c' M } for H2 (K3, Z) is defined to have intersection matrix 

~U 



h'MN — #( C M) c 'n) 



u 



U 



-Eh 



(B.7) 
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where E% represents the Cartan matrix of that group, 

-2 1 
1 -2 1 
1 -2 1 

— Est = 

1-211 
1 -2 

1 -2 1 
1 -2 

The change of basis relating the {c/} to the {c' M } is as follows: 

c' M = M M T ci, 

where 



2M M ! = 
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(B.10) 

(Thus one has h' MN = Mm 1 M^ J hu .) In deriving this change of basis we made use of results 
in |l3|| . To our knowledge its explicit form is new. 



(B.8) 
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